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We report observations of CP violation in the decays B° — > K + ty and B° — > 7r + 7r in a sample 
of 383 million r(4S) -> BB events. We find 4372 ± 82 B° -> K + -k~ decays and measure the 
direct CP-violating charge asymmetry Aktt = —0.107 ±0.018 (stat)^'^ (syst), which excludes the 
CP-conserving hypothesis with a significance of 5.5 standard deviations. In the same sample we 
find 1139 ± 49 B° — > 7r + 7r~ decays and measure the CP-violating asymmetries = —0.60 ± 
0.11 (stat) ±0.03 (syst) and C™ = -0.21 ±0.09 (stat) ±0.02 (syst). CP conservation in B° -» tt+tt" 
{Stttt = = 0) is excluded at a confidence level 1 — C.L. = 8 x 10~ 8 , corresponding to 5.4 standard 
deviations. 

PACS numbers: 13.25.Hw, 11.30.Er, 12.15.Hh 



The prediction of large CP-violating effects in the B- 
meson system [l[ has been confirmed in recent years by 
the BABAR and Belle collaborations, both in the inter- 
ference of B decays to charmonium final states with 
and without B°-B° mixing [2|, and directly in the 
interference between the decay amplitudes in B —>■ 
All measurements of CP violation to 
date are in agreement with indirect predictions from 
global standard-model (SM) fits [6| based on measure- 
ments of the magnitudes of the elements of the Cabibbo- 
Kobayashi-Maskawa (CKM) quark-mixing matrix [7J 
and place important constraints [8[ on the flavor struc- 
ture of SM extensions. 

The proper-time evolution of the asymmetry between 
B° and B° decays to ir + ir~ is characterized by sine 
and cosine terms with amplitudes Sim, which arises 
from interference between decays with or without B°- 
B° mixing, and C^, which is due to interference be- 
tween the b — > u "tree" and the higher-order b — > d 
"penguin" decay amplitudes. Similarly, the direct-CP- 
violating asymmetry Ak-k between the B° — > K~ir + and 
B° — > K + ir~ decay rates arises from interference between 
b — ► u tree and b — > s penguin amplitudes. Negligi- 
ble contributions to these asymmetry parameters would 
also enter from CP violation purely in B°-B° mixing, 
which has been determined to be very small The 
quantity sin2a c ff = SV^/ yl — C^ can be related to 
a = arg [— V td V t ^/V ud T/ u * b ] through a model-independent 
analysis that uses the isospin-related decays B ± — > 



and B 



o 



Contributions from new particles 



could affect the asymmetries in these modes primarily 
through additional penguin B-decay amplitudes. 

Previous evidence of direct CP violation in B° — > 
K + ir~ has been reported by BaBAR Q and Belle ad- 
ditional measurements of Aktt have also been reported 
by the CDF [ll| and CLEO [H Collaborations. The 
Belle Collaboration recently reported [l3| an observa- 
tion of both time-dependent and direct CP violation in 
B° — > 7r + 7r~ decays using a sample of 535 x 10 6 BB 
pairs, while our previous measurement 14| on a sample 
of 227 x 10 6 BB pairs was statistically consistent with no 
CP violation. In this Letter, we present measurements 
of Ak-w, Sktt, and C^ in a sample of 383 x 10 6 BB pairs 
using an improved analysis technique with significantly 
increased sensitivity compared to our previous measure- 
ments. 



In the BaBAR detector [15J, charged particles are de- 
tected and their momenta measured by a combination of 
a five-layer silicon vertex tracker (SVT) and a 40-layer 
drift chamber (DCH) that covers 92% of the solid angle 
in the T(4S I ) center-of-mass (cm.) frame, both operat- 
ing in a 1.5-T solenoidal magnetic field. Discrimination 
between charged pions, kaons, and protons is provided 
by a combination of an internally reflecting ring-imaging 
Cherenkov detector (DIRC), which covers 84% of the cm. 
solid angle in the central region of the BABAR detector 
and has a 91% reconstruction efficiency for pions and 
kaons with momenta above 1.5 GeV/c, and the ioniza- 
tion (dE/dx) measurements in the DCH. Electrons are 
explicitly removed based on a comparison of the track 
momentum and the associated energy deposition in a 
CsI(Tl) electromagnetic calorimeter, and with additional 
information from dE/dx and DIRC Cherenkov angle (8c) 
measurements. 

The analysis method retains many features of our 
previous B° — > K + tt~ and B° — > tt + it~ CP-violation 
measurements 0, [l4| • We reconstruct candidate decays 
B rC c — > h + h~ (h ± = 7r ± , K ± ) from pairs of oppositely 
charged tracks in the polar-angle range 0.35 < 8\ a b < 
2.40 that are consistent with originating from a common 
decay point. The remaining particles are examined to 
infer (flavor tag) whether the other B meson in the event 
(•Stag) decayed as a B° or B°. We perform an unbinned 
extended maximum-likelihood (M.L.) fit simultaneously 
for the CP-violating asymmetries and the signal and 
background yields and parameters. The fit uses particle- 
identification, kinematic, event-shape, B tag flavor, and 
At information, where At is the difference between the 
B roc and B tag decay times. The yields for the Kit final 
state are parametrized as n K ± nT = n K „ (1 qp -^kZ) 
and the decay-rate distribution / + (/_) for B roc — ► 7r + 7r~ 
and B t ag = B° (B°) is given by 

e -|At|/r 

/±(At) = — [1 ± sin(Am d At) 

4r 

T C 7r7r cos(Am (J At)], (1) 

where r is the neutral B lifetime and Am^ is the B°-B° 
mixing frequency, both fixed to their world averages [§]. 

The most significant improvement in sensitivity com- 
pared to our previous analysis comes from a 35% increase 
in the B rec reconstruction efficiency that results from us- 
ing dE/dx as a discriminating variable in the M.L. fit 
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for the first time. The dE/dx measurements are used 
both to complement the discriminating power of 6q for 
charged particles within the DIRC acceptance and as 
a standalone means of particle identification for tracks 
that have no DIRC information and were not included 
in our previous measurements. The dE/dx calibration 
takes into account variations in the mean value and res- 
olution of dE/dx with respect to changes in the DCH 
running conditions over time and each track's charge, po- 
lar and azimuthal angles, and number of ionization sam- 
ples. The calibration is performed with large (> 10 6 ) 
high-purity samples of protons from A — > p7r~, pions and 
kaons from D* + — > D°tt + (D° — > K~tt + ), and additional 
samples of pions from r~ — > 7r~7r + 7r~^ r decays and from 
K® — > 7r + 7r~ decays that occur in the vicinity of the in- 
teraction region. 

We require at least one of the tracks in the B lec de- 
cay candidate to have 9c measured with at least six sig- 
nal photons; for such tracks, the value of 9c must agree 
within 4 standard deviations (a) with either the pion or 
kaon hypothesis. Thus, protons with six or more sig- 
nal photons are removed while proton-pion and proton- 
kaon combinations are possible for background candi- 
dates where one of the tracks has no usable 9c mea- 
surement. We construct 9c probability-density functions 
(PDFs) for the pion and kaon hypotheses, and dE/dx 
PDFs for the pion, kaon, and proton hypotheses, sepa- 
rately for each charge. The K-tt separations provided 
by 9c and dE/dx are complementary: for On, it varies 
from 2.5a at 4.5GeV/c to 13cr at 1.5GeV/c ff, while for 
dE/dx it varies from less than l.Ocr at 1.5 GeV/c to 1.9a 
at 4.5GeV/c (Fig.H). 

Each B candidate is characterized by the energy dif- 
ference AE = (qr ■ qs/y/s) — yfs/2, which also provides 
additional discriminating power between the four pos- 
sible final states (ir + ir~, K + ir~, K~tt + , K + K~), and 
the beam-energy-substituted mass Toes = [( s /2 + fir ■ 
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[151 ]. Here, qr and qs are the four- 
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FIG. 1: The average difference between the expected val- 
ues of DIRC 8c and DCH dE/dx for pions and kaons at 
0.35 < #i a b < 2.40, divided by the uncertainty, as a function 
of laboratory momentum in B a — » h + h~ decays in BABAR. 



momenta of the T(AS) and the B candidate, s = {qr) 2 
is the square of the cm. energy, fir and pb are the 
laboratory three-momenta of the T(4S) and the B, and 
Er = q T is the laboratory energy of the T(45). For signal 
events, the toes and AE PDFs are Gaussian functions 
with widths of 2.6MeV/c 2 and 29MeV, respectively. For 
the background, tors is parametrized with an empirical 
threshold function [l6| and AE is parametrized with a 
second-order polynomial. Wc require 5.2 < toes < 5.3 
GeV/c 2 and \AE\ < 0.150 GeV. 

The background arises predominantly from random 
combinations of tracks in e + e _ — > qq (q — u, d, s, c) and 
t + t~ jetlike continuum events. We define the angle 9s in 
the cm. frame between the sphericity axes 17[ of the B 
candidate and of all remaining charged and neutral par- 
ticles in the event. For background events, |cos6*s| peaks 
sharply near 1, while for B decays the distribution is 
nearly flat. We require |cos6>s| < 0.9, which removes ap- 
proximately 64% of uu, dd, and ss, 52% of cc, and 84% of 
t + t~ background. Contamination from e + e~ — * t + t~ 
production is reduced to 2% of the total background by 
requiring the ratio of the second to zeroth Fox- Wolfram 
moments 18j to be less than 0.7, which has a negligible 
effect on the signal efficiency. The overall gain in signal 
reconstruction efficiency is 52% compared to our previous 
analysis. Additional continuum-background suppression 
in the fit is accomplished by the Fisher discriminant T 
described in Ref. [l9(. We have studied the backgrounds 
from higher-multiplicity B decays and find them to be 
negligible, particularly due to their good separation from 
signal in AE. 

The Stag flavor is determined with a neural-net algo- 
rithm [2(| that assigns the event to one of seven mu- 
tually exclusive tagging categories. The figure of merit 
for the tagging quality, measured in a data sample -Bflav 
of fully reconstructed B° decays to Z?W~(7r + , p+, a+) or 
J/ipK*°, is the effective efficiency Q = J2k e k(^ — ^w k ) 2 — 
0.305 ± 0.003, where e k and Wk are the efficiencies and 
mistag probabilities for events in tagging category k. 
Separate values of e k and Wk for each background cat- 
egory are determined in the M.L. fit. 

The time difference At = Az/fijc, where (3j ~ 0.56 
is the known boost of the T(4S f ), is obtained by measur- 
ing the distance Az along the beam (z) axis between the 
B rcc and B tag decay vertices. We require |Ai| < 20 ps 
and (JAt < 2.5 ps, where (JAt is the At uncertainty esti- 
mated separately for each event. The resolution function 
for signal candidates is a sum of three Gaussians [2(| 
with parameters determined from a fit to the full -Bflav 
sample. The background At distribution, common to all 
tagging categories, is modeled as a sum of three Gaussian 
functions with parameters determined in the final fit. 

The likelihood for candidate j tagged in category k is 
obtained by summing the product of event yield r^, tag- 
ging efficiency e^.fc, and probability Vi.k over all possible 
signal and background hypotheses i. We treat separately 
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the cases where both or only one track has a 6c mea- 
surement. The extended likelihood function for tagging 
category k is 



Ck = cxp - ^2 n id,k J \ 



y^ j n l e it k'PiM{xj;d l ) 



The probabilities V^k are evaluated 
PDFs for each of the independent 
{toes, AE, T, dE/dx, 6c, At}, with parameters on 
use separate 6q and dE/dx PDFs for positively 



(2) 

as a product of 
variables Xj = 
We 
and 



negatively-charged tracks. The At PDF for signal 7r + 7r~ 
decays is given by Eq. ([T|) modified to include the mistag 
probabilities for each tagging category and convolved 
with the signal resolution function. The At PDFs for sig- 
nal Ktt and background Ktt, irp, and Kp combinations 
take into account the correlation between the charge of 
the kaon or proton and the -Btag flavor; for signal Ktt, 
B°-B° mixing is also taken into account. The total like- 
lihood C is the product of likelihoods for each tagging 
category and has 117 free parameters. 

Fitting the final sample of 309540 events, we find 
= 1139 ± 49, n K „ = 4372 ± 82, n KK = 10 ± 17, 
where all errors are statistical only, and measure the fol- 
lowing asymmetries: 

Aktt = -0.107 ± 0.018 (stat)+°;°°^ (syst), 
S„„ = -0.60 ± 0.11 (stat)± 0.03 (syst), 
CW = -0.21 ± 0.09 (stat) ±0.02 (syst). 

Here Ak-k is the fitted value of the K^tt^ event-yield 
asymmetry A T g™ shifted by +0.005±j^ to account for 
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a bias that arises from the difference between the cross 
sections of K + and K~ hadronic interactions within the 
BABAR detector. We determine this bias from a detailed 
Monte Carlo simulation based on GEANT4 [2l[ version 
7.1; it is independently verified with a calculation based 
on the known material composition of the BABAR de- 
tector |15| and the cross sections and material proper- 



ties tabulated in Ref. @. The corrected K^tt^ event 
yield asymmetry in the background, where no observ- 
able CP violation is expected, is consistent with zero: 
-0.006 ± 0.004 (stat)t°;°°s ( syst ). 

A contour plot of the C ffff ) confidence levels is 

shown in Fig. [2j The correlation between and C vv 
is —0.07. Performing a fit that excludes At and using an 
event- weighting technique 22j, in Fig. [3] we show the dis- 
tributions of At for signal 7r + 7r~ events with -Btag tagged 
as B° or B°, and the asymmetry as a function of At, 
overlaid with the PDF curves that represent the result of 
the full fit. 

To validate our results, we perform a number of consis- 
tency checks and systematic-error studies similar to those 
reported in Refs. Q and 14 1. For Ak-k, the dominant 



source of systematic uncertainty is the bias due to kaon 
hadronic interactions. We find that the systematic errors 
due to potentially imperfect understanding of the DIRC 
and DCH particle-identification performance are small 
iovA K * (0.002), (0.007), and (0.006). The dom- 
inant sources of systematic uncertainty on S^k are the 
signal At model (0.020) and flavor-tagging parameters 




FIG. 2: SVtt and C ww : the central values, errors, and 
confidence-level (C.L.) contours for 1 — C.L. = 0.317 (lcr), 
4.55xl0~ 2 (2a), 2.70xl0 -3 (3a), 6.33xl0~ 5 (4cr), 5.73xl(T 7 
(5<j), and 1.97 x 10~ 9 (6<r), calculated from the square root 
of the change in the value of — 21n£ compared with its value 
at the minimum. The systematic errors are included. 



FIG. 3: The background-subtracted distributions of the 
decay-time difference At in signal B — > 7v + n~ events. The 
points with errors show the events where -B tag is identi- 
fied as (a) B° or (b) B°. The asymmetry, defined as 
(n B o — ngo) / (n B o + rigo), for signal events in each At bin, 
is shown in (c). The solid curves are the projection of the fit. 
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(0.015), while for C n7r the dominant uncertainties arise 
from tagging (0.014) and the potential effect Q of dou- 
bly CKM-suppressed decays of the B tag meson (0.016). 
As a final cross-check, we perform a fit allowing the mix- 
ing frequency and lifetime to vary simultaneously with 
SW and CU. Wc find Am d = 0.506 ± O.O^ps" 1 and 
t b o = 1.523 ± 0.026 ps, where the errors are statistical 
only, consistent with the world-average values, and the 
resulting shifts in the CP parameters arc negligible. The 
total systematic uncertainties are calculated by summing 
all individual contributions in quadrature. 

In summary, we observe direct CP violation in the de- 
cay B° — > K + 7r~ with a statistical significance of 5.5a 
and CP violation in the time distribution of B° — > ir + ir~ 
decays with a significance of 5.4er. We also determine 
that the mixing-induced CP-violating asymmetry S vn is 
nonzero with a significance of 5.1er or greater for any 
value of Ctt^. All results are consistent with, and super- 
sede, our previously published measurements 0, 
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